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Abstract. Building upon a previous analysis of P v wind lines in LMC O 
stars, we analyze the P v wind lines in a sample of Galactic O stars which 
have empirical mass loss rates determined from either their radio fluxes or Ha 
profiles. Since the wind analysis provides a measure of Mq where q is the 
ionization fraction of the ion, we determine q(P v) observationally. In spite of 
model predictions that q ~ 1 for mid-0 stars, we find q(P v) < 0.15 throughout 
the O stars. We discuss the origin of this discrepancy. 



1. Measures of M 

Three approaches are normally used to determine stellar mass loss rates, Ms. 
All assume that the wind is homogeneous and spherically symmetric (SS) with a 
single, monotonic velocity law, and all should agree. The three approaches are: 

1. Continuum excess from free- free emission. This samples the outer wind 
(the exact radius depends on wavelength), where it becomes optically thick to 
free-free emission. It is only detectable for massive winds in nearby stars. The 
radio wavelengths are considered "cleanest", because, in contrast to the IR/FIR, 
massive winds become optically thick at large radii (> 10i?*), where v = Voo, a 
constant. This makes p w i n d ~ M /(r 2 Voo), independent of v(r). Furthermore, no 
photospheric correction is needed. However, the radio flux can be non-thermal, 
so observations at multiple wavelengths are required to determine the spectral 
index of the emission. 

2. Ha emission. This samples the inner wind and is easily observed. For massive 
winds, Ha emission is related to M. However, the exact form of the observed 
Ha profile depends upon the N = 3 departure coefficient for H in the wind 
and this, in turn, depends upon: the photospheric radiation field; the diffuse 
radiation field of the wind; and the wind velocity law in the acceleration region. 
The shape of the "photospheric" Ha profile is also required, and the observed 
W\{Ha) can be strongly variable. Nevertheless, relatively sophisticated models 
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for Ha formation exist (e.g., Repolust et al. 2004), and can provide reasonable 
agreement between available radio and Ha Ms when Ha emission is strong. 
3. UV resonance lines. These sample the entire wind. Their shapes are deter- 
mined by the radial optical depth of the wind, T ra d ~ MqiAE, where Ae and qi 
are the abundance of element E and its ionization fraction for stage i. However 
observations of a dominant ion (qi ~ 1), of known abundance are required to es- 
timate M directly, and the wind lines of abundant, dominant ions are saturated 
in winds massive enough to be detected in the radio or to have reliable Ha Ms. 

2. FUSE and P V 

FUSE gives access to P v AA1118, 1128. P v is a surrogate for C iv (Massa et 
al. 2003) and qi ~ 1 is expected for both ions in mid-0 star winds. Unlike Ac, 
Ap ~ Const over the life of an O star. Furthermore, for scaled solar abundances, 
Trad(C iv)/r m rf(P v) = 661, so to detect P v, T ra d(C iv) > 50 and saturated - 
as is the case for stars with detectable radio fluxes and strong Ha emission. 

P V in LMC O Stars: The first large scale FUSE study of P v was by Massa 
et al. (2003). They performed SEI (Lamers et al. 1987) fits to P v wind lines, 
determining T ra d for 25 LMC O stars. They then used Ms predicted by the Vink 
et al. (2001) theory and Ap = 0.50 x solar to find that q{P V) peaked between 
45-50 kK, as expected, but with a peak value < 0.15 (see Fig. 1). This result 
implies a factor of 7 or more discrepancy between the expected and observed 
Ms. There are three possible explanations for this result: 

1. The LMC P abundance scales differently from other elements, 

2. The theoretical Ms are incorrect for the LMC, or 

3. The winds are not homogeneous and SS, but strongly clumped or structured. 

P V in Galactic O Stars (preliminary results): The Galactic Ap is well 
determined (e.g., Catanzaro et al. 2003), so abundance is not an issue for the 
Galaxy. We are currently analyzing P v in Copernicus, Orfeus and FUSE data 
for stars with radio and/or Ha M estimates. This eliminates the need for model 
Ms to derive q for P v, and determines how well different mass loss rate indi- 
cators agree. So far, we have analyzed 30 stars (11 more will be observed by 
FUSE), and the results are shown in Figure 2. These results are effectively 
identical to those from the LMC (Fig. 1), implying the same conclusions. 

3. Conclusions 

• An erroneous P abundance is not the cause of the small LMC q(P V)s. 

• Strong clumping/porosity must be the root of the problem. 

• Large scale clumping can also strongly affect the radio and Ha Ms, so 
these should be re-evaluated. 

• The good news: the different measures of M are sensitive to different as- 
pects of the wind flow, so bringing them all into agreement (together with 
the X-ray and O vi wind line observations) will provide powerful con- 
straints on models of how the winds are structured. 
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Figure 1. Mean P v ionization fractions as a function of temperature 
for LMC O stars analyzed by Massa et al. (2003). Open, half filled and 
filled symbols denote stellar luminosities in the ranges logL/L© > 6.0, 
6.0 > logi/L > 5.6 and 5.6 > logL/L Q , respectively. 
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Figure 2. Mean P v ionization fractions versus temperature for Galactic O 
stars with radio or Ha Ms. Filled symbols: radio Ms, open circles: Ha Ms 
from Repolust et al. (2004), open triangles: Ha Ms from Lamers & Leitherer 
(1993). 
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